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We consider a Bohmian approach to the Wheeler-DeWitt quantization of the Friedmann-Lemaitre- 
Robertson-Walker model and investigate the question whether or not there are singularities, in the 
sense that the universe reaches zero volume. We hnd that for generic wave functions (i.e., non- 
classical wave functions), there is a non-zero probability for a trajectory to be non-singular. This 
should be contrasted to the consistent histories approach for which it was recently shown by Craig 
and Singh that there is always a singularity. This result illustrates that the question of singularities 
depends much on which version of quantum theory one adopts. This was already pointed out 
by Pinto-Neto et ai, albeit with a different Bohmian approach. Our current Bohmian approach 
agrees with the consistent histories approach by Craig and Singh for single-time histories, unlike the 
one studied earlier by Pinto-Neto et al. Although the trajectories are usually different in the two 
Bohmian approach, their qualitative behavior is the same for generic wave functions. 

PACS numbers: 98.80.Cq, 04.60.Ds, 98.80.Qc 


I. INTRODUCTION 

Recently, some of us investigated the issue of singular¬ 
ities in the Wheeler-DeWitt approach to quantum cos¬ 
mology and showed that the answer depends strongly on 
the version of quantum mechanics that one considers [T] . 
Two versions were compared: the consistent (or deco- 
herent) histories and the Bohmian approach (also called 
the deBroglie-Bohm or pilot-wave approach). Both ap¬ 
proaches have a way of introducing possible histories of 
the universe. In the Bohmian approach, there is an ac¬ 
tual configuration whose time evolution is determined by 
the wave function. The possible paths of the configura¬ 
tion are the histories. On the other hand, in the consis¬ 
tent histories approach the (coarse-grained) histories are 
sequences of propositions at different moments in time. 
(For example, a possible proposition could be that the 
system occupies a certain region in space). 

A simple quantum cosmological model was studied, 
namely the quantized Friedmann-Lemaitre-Robertson- 
Walker model, with a homogeneous scalar field. For 
the consistent histories approach, Craig and Singh [ai 
showed that the possible histories are always singular ir¬ 
respective of the given wave function. Either they start 
from a singularity or they end in a singularity. On the 
other hand, according to the Bohmian approach of Pinto- 
Neto et al. [T], there are wave functions for which there 
are also histories that display a bounce (which corre¬ 
sponds to a universe that contracts until it reaches a 
minimum volume and then starts expanding) and hence 
do not have a singularity. 
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However, while the dynamics of the Bohmian model is 
very natural, there is no natural probability distribution 
on the set of histories. As such, there is no immediate way 
of making probabilistic statements, like about the prob¬ 
ability for a history to have a singularity for a particular 
wave function. On the other hand, Craig and Singh laE] 
showed that probabilistic statements can be made in the 
consistent histories approach. This might seem puzzling 
since, at least in the context of non-relativistic systems, 
it is often claimed that consistent histories and Bohmian 
mechanics yield the same predictions for outcomes of 
measurements IS 0. The reason for this mismatch is 
that in the consistent histories approach the scalar field 
is treated (at least formally) as a time variable, whereas 
it is not in the Bohmian approach. The goal of this pa¬ 
per is to explore another Bohmian dynamics which is also 
based on treating the scalar field as a time variable. In 
this case, there is a natural probability distribution for 
the scale factor which agrees with the one of the con¬ 
sistent histories approach when single-time histories are 
considered. Hence, this allows for a better comparison 
between the Bohmian approach and the consistent histo¬ 
ries approach concerning the question of singularities. In 
particular, given a wave function, we can now calculate 
the probability for a trajectory to run into a singularity. 
We will find that although the trajectories are different 
from the ones in the model proposed in [T], they are of¬ 
ten qualitatively the same. In particular, we will find that 
some wave functions allow for trajectories that are never 
singular and instead display a bounce. More precisely, 
we will see that for a non-classical wave function there is 
always a non-zero probability for a bounce (while trajec¬ 
tories corresponding to a classical wave function always 
have a singularity). This should be contrasted with the 
fact that in the consistent histories approach the proba¬ 
bility for a singularity is always one, for any wave func¬ 
tion. 
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The outline of the paper is as follows. We start 
by recalling the Wheeler-DeWitt quantization of the 
Friedmann-Lemaitre-Robertson-Walker model and the 
Bohmian approach of [T]. In section III we will discuss 
possible choices of Hilbert space that were considered in 
the consistent histories approach. In section |IV[ we will 
develop and discuss the alternative Bohmian dynamics. 
In section|Vj we calculate the probability for a trajectory 
to be singular. Finally, in section |VI[ we consider the 
probabilities for singularities in the consistent histories 
approach, extending the work of Craig and Singh from 
two-time histories to many-time histories, and compare 
these to those predicted by the Bohmian approach. 


II. WHEELER-DEWITT QUANTIZATION AND 
BOHMIAN MECHANICS 


A classical flat Friedmann-Lemaitre-Robertson-Walker 
space-time is described by a metric 

ds^ = , (I) 

where N is the lapse function, a = e“ is the scale fac¬ 
tor, and 'jij is the flat spatial 3-metric. Considering the 
matter content described by a free massless homogeneous 
scalar held the Lagrangian^ of the system can be 
written as [5] 


where a dot means derivative with respect to coordinate 
time t. The corresponding equations of motion lead to 


= ±- 


N 

,3a 


C, 


a = 


p3q 


( 3 ) 


where c is an integration constant. N remains an arbi¬ 
trary function of time. This implies that the dynamics 
is time reparameterization invariant. Different choices of 
N merely correspond to different choices of the parame¬ 
terization of the paths in (i^, a) space. In the case c = 0, 
the universe is static. For c ^ 0, we have 


a = dz(j) -I- c, 


( 4 ) 


with c another integration constant. 

The universe reaches the singularity (i.e., zero volume) 
when a = e“ = 0. For a non-static universe the singular¬ 
ity is reached for either (j) —> —oo or 0 —>■ oo. So the uni¬ 
verse either starts with a big bang or ends in a big crunch. 
In terms of proper time r, which is defined by dr = Ndt, 


^ For simplicity, we have dropped the constant factor ATrGfS in 
front of a^. This factor could be removed by a suitable rescaling 
of the scalar field. 


integration of ([^ yields a = e“ = [ 3 (ct -|-where 
c is an integration constant, so that a = 0 for t = —c/c 
(and there is a big bang if c > o and a big crunch if c < 0). 
This means that the universe reaches the singularity in 
finite proper time. 

Canonical quantization of this theory leads to the 
Wheeler-DeWitt equation 

dl^-dli; = 0. ( 5 ) 

In the corresponding Bohmian theory there is an ac¬ 
tual scalar field (j) and scale factor a = e“, which satisfy 

« = ( 6 ) 

where il) = I'i/'je''^. The function N is again the lapse 
function, which is arbitrary, and, as in the classical case, 
implies that the dynamics is time reparameterization in¬ 
variant.^ 

The Wheeler-De Witt equation implies 

{d^Sf-{d^Sf + Q = 0, (7) 


{\ijj\‘^d^S) - da {\tjj\'^daS) = 0, ( 8 ) 


where 


0 = + 


(9) 


is the quantum potential. If Q = 0, then Q implies that 
= {daS)"^. In addition, daS is then conserved 
along a Bohmian trajectory. Hence, in this case, the 
Bohmian motion is reduced to the classical motion given 
by 

Equation Q is a continuity equation and implies that 
the Bohmian dynamics preserves However, since 

\ip\'^ is not normalizable it can not be straightforwardly 
used to make statistical predictions. In this paper, we 
will consider an alternative Bohmian dynamics, which 
allows for immediate statistical predictions and which 
allows for a direct comparison with consistent histories 
approaches discussed in mmi]. 


^ Note that the time reparameterization invariance is a special fea¬ 
ture of Bohmian approaches to mini-superspace models [3 ig. 
For the usual formulation of Bohmian dynamics for the full 
Wheeler-DeWitt theory of quantum gravity, a particular space¬ 
like foliation of space-time or, equivalently, a particular choice 
of “initial” space-like hypersurface and lapse function, needs to 
be introduced. Different foliations (or lapse functions) yield dif¬ 
ferent Bohmian theories [9ini|. So, in this case, the dynamics 
is not invariant under space-time diffeomorphisms. (Yet, while 
the usual Bohmian formulation is not diffeomorphism invariant, 
this may perhaps be achieved with alternative approaches. For 
a discussion of analogous issues concerning special relativity in 
Bohmian mechanics, see m) 
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which satisfy^ 

=H±'il)± = T\/-^V'± ■ 


(15) 



a 



FIG. 1: The wave functions tpR and ipL respectively move to 
the right and to the left without change in shape. 


III. WHEELER-DEWITT EQUATION AND 
HILBERT SPACES 


The left- and right-moving components of ^p± will be de¬ 
noted by '4’l,± and 4’R,±- If one uses the scalar held (j) as 
a time variable, equation (151 represents a Schrddinger- 
like equation and given the initial wave function at 
time 00 j tlio (formal) solution is given by 


0±(a, 0) = e^'V^ ‘^“^0^(a). (16) 

The ftrst choice of Hilbert space is "Hi = L^(M,C^) 
m [31 Hi]. We have that %i — 'Hi,+ 0 "Hy-, with 
'Hi^± = L^(M,C) the positive and negative frequency 
Hilbert spaces. We write 'F(a) = (0+(a),0_(a)) with 
'F G "Hi and 0± G The inner product is given by 


So far we have not introduced a Hilbert space for the 
Wheeler-DeWitt equation ([^. We will discuss two pos¬ 
sible choices of Hilbert space that have appeared in the 
literature jaiidiHiii- Both are motivated by consid¬ 
ering 0 (at least formally) as a time variable. However, 
as we shall see, for the purpose of assigning histories to 
the scale factor, either in the consistent histories or in 
the Bohmian framework, both choices can be considered 
equivalent, assuming a suitable choice of “observable” for 
the scale factor in each case. 

First, consider a general solution to the Wheeler- 
DeWitt equation ([^, which is of the form 

0(a, 0) = 0L(a + 0) + 0fi(a - 0), (10) 

where the indices L and i? denote respectively “left- 
moving” and “right-moving”. This terminology stems 
from the fact that if 0 is regarded as time, then ipL, as 
a function of a, moves to the “left” over time, without 
change in shape, whereas tpR moves to the right, see fig¬ 
ure In terms of the Fourier transform, we have 

0(0,0)=^/ dA 0L(fc)e‘'=(“+« 

V 27r J—QQ 

1 

+ — / dk 0fl(A)e‘'=(“-^). ( 11 ) 

V 27^ J — oo 

The general solution can also be decomposed into a pos¬ 
itive and negative frequency part 

0 = 0++0_, (12) 

where 

1 

0+ / dk 0i(fc)e''=(“+'^) 

v27r Jo 

+ ^ r dk 0i^(fc)e‘'=(“-^), (13) 

V 27r J —oo 

1 

0_ =-= / dk 0i(A;)e‘'=(“+'^) 

V J — QQ 

+ ^ / dk 0fl(fc)e‘'=(“-« , (14) 

v27r Jo 


(vI/|ci>)i = y'da(0;0+ + 0*0_) . (17) 

The dynamics for the positive and negative frequency 


states is given by (15). Since the Hamiltonians H± are 


Hermitian operators on this dynamics preserves 

the inner product. 

The natural observable corresponding to the scale fac¬ 
tor is given by the projection-valued measure (PVM) 


Pi(da) = Pi,+ (da) 0 Pi,_(da), (18) 

where 

Pi,±{da) = |a)(a|da, (19) 

with jd) = 6{a — a). (Since f Pi^±(da) = 1±, with 
the identity operator on 'Hi^±, it immediately follows that 
/ Pi (da) = 1, the identity operator on Pi.) 

For a state 4' = (0+,0_), the corresponding density 
for a is 


p(a,0) = |0+(a,0)p 0 |0_(a,0)p , (20) 

i.e., (4'|Pi(da)|'I')i = p{a,(j))da. (The Hilbert space is 
the same as that of the Pauli theory for a non-relativistic 
spin-1/2 particle, with Pi (da) being the analogue of the 
position PVM in the Pauli theory.) 

The second choice of Hilbert space P2 is again the 
direct sum of positive and negative frequency Hilbert 


® For a wave function the 

action of the Hamiltonian is defined as 9^ 1 / 1 ( 0 ) = 

^ Actually, in [lEllTl] only the positive frequency sector of "Hi is 
considered, by appeal to superselection. However, the Bohmian 
trajectories will depend on whether the wave function is a super¬ 
position of both frequencies or not. (Following the terminology 
of |17| . one could say that frequency is weakly superselected, but 
not strongly.) Therefore, a general analysis should include both 
frequency sectors. 
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spaces, i.e., 'H 2 = 'H 2 + ©"^ 2 -, with 'H 2 + = 'H 2 - 
m [16]. States in 'H 2 are maps from M to and 
states in 'H 2 ,± are maps from M to C. We write ^(a) = 
('0+(a),'0_(a)) with ^ G'H 2 and 'ip± G 7^2,± (the bar is 
used in order to distinguish the states from those in "Hi). 
The inner product on H 2 ,± is given by 

(^l<^)2.± = 2 J da-ilj*{TH±)^ = 2 J da^* , 

( 21 ) 

so that the inner product for "^2 reads® 

(§|$)2=2 I da(^;^/^<^+ + C^/=^^-) • (22) 


The dynamics for the positive and negative frequency 
states is given again by (151. Since the Hamiltonians 


H± are Hermitian operators on 712,±, the dynamics also 
preserves this inner product. 

For the norm of a wave function 'h, the integrand on 
the right hand side of (22) reads 


2 (iA+\/^V^+• 


(23) 


One might be tempted to regard this quantity as the 
density of a at a given time (p for the state However, it 
is not always positive dehnite, even if we restrict ourselves 
to positive frequencies [Ulllg]. 

One possibility to obtain a positive density is by using 
the Newton-Wigner states (TU |5D| 


\a,NW) 



^ik(a-a) 

Vm 


(24) 


These states are normalized according to 
{a,NW\a,NW )2 = 5{a — a) and form a basis of 
7^2,±- The PVM then reads 


P2{da) = P2^+{da) ® P2-{da), (25) 


where now 


P2,±{da) = \a,NW){a,NW\da. 


(26) 


Note that, as before, / P 2 {da) = 1, being understood 
that the action of these operators on the states in the 
Hilbert space 772 is through the inner product defined in 
Eq. (|22l). 


® Perhaps a more familiar way of writing the inner products | |21| l 
and | |22| l is in terms of solutions to the wave equations | |15| l 1141 . 
The inner product ( |21[ | then corresponds to the Klein-Gordon 
inner product {'P\4>)kg = for two pos¬ 

itive frequency solutions to ( |l5| |. The inner product | |22| l then 
corresponds to (4'|4>)2 = {b-tlJb+)ifG ~ iP—lP—)KG- Note that 
this inner product is positive definite, unlike the usual definition 
for the inner product of two solutions to the Klein-Gordon equa¬ 
tion, which does not contain the relative minus sign between the 
positive and negative frequency part. 


In the consistent histories approach, (coarse-grained) 
histories for a can be introduced, by choosing a Hilbert 
space, a PVM (or positive operator-valued measure 
(POVM)) and a Hamiltonian. The results will be iden¬ 
tical when choosing either (77i,Pi,.f7) or {'H 2 , P 2 , H), 


where H is the Hamiltonian operator determined by (15) 
(the former is considered in [31|3], albeit just for positive 
frequencies, the latter in [H]).® The reason is that there 


is a unitary mapping 7^2 —t 77i, which maps P 2 -G Pi 
and leaves the Hamiltonian invariant. This mapping 
is defined by 'ip±{a) —>■ i/’±(a) = {a, NW\p)±) 2 ,± (or 
\a,NW) —>■ |a) in terms of basis states, or 'ip±{k) —>■ 
'tp±{k) = •^2|fc|'i/j±(/c) in terms of Fourier components). 
This map is unitary since (4'|<i))i = (^|4>)2 (which follows 
from the fact that f P 2 (da) = 1, the identity operator 
on 7 ^ 2 ). Since it takes P 2 —>■ Pi, we have 


('I'|P2(rfa)|'k)2 = ('k|-Pi(rfa)|'k)i = p{a,(j))da . (27) 


As such, we clearly get the same result for probabilities 
of histories with (77i,Pi) or (7^2, P 2 )- 

We can also develop a Bohmian approach given a 
Hilbert space, a POVM and a Hamiltonian. For a 
given wave function, the distribution corresponding to 
the PVM will be the Bohmian equilibrium distribution. 
The dynamics can then be defined in such a way that this 
distribution is preserved [221 [23] • In the next section, 
we will consider a dynamics that preserves the density 
(20). So also in the Bohmian approach, the two choices 
{Til, Pi, H) and {712, P 2 , H) yield an identical dynamics 
for the scale factor. 


IV. ALTERNATIVE BOHMIAN DYNAMICS 


The density p{a,(p), given in (20), is not preserved by 
the Bohmian dynamics (§,i.e., if the density of a is given 
by p{a, po) at a certain time po, then the Bohmian dy¬ 
namics in general leads to a density different from p{a,p) 
at other times p [Tj. However, one can postulate a differ¬ 
ent dynamics that preserves p [23U2H] . According to this 
dynamics, possible trajectories a{p) satisfy 


ra(4>) 


dap{a, p) = 


l-a{0) 


dap{a, 0), 


(28) 


or, equivalently, 

dp p{a,p)J_^ 


dad^p{a, p). (29) 


Since p satisfies the continuity equation 

d^p + da {V2P) = 0 , 


(30) 


See [21| for similar considerations on equivalence of quantum 
theories. 
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it is preserved by the Bohmian dynamics. 

Let us compare this dynamics to the one formulated 
in section [T^ First, consider a real wave function ip, i.e. 
Im Ip = 0. According to the dynamics of section |nj 
the possible solutions are a{t) = a(0), p{t) = ^(0), so 
that the universe is static, with a constant scale factor 
and constant scalar field. Since these solutions can not 
be expressed as a{p) they must be different from the 
trajectories given by (28). 


Now assume that the trajectories for the first Bohmian 
dynamics can be expressed as a{(j)) (at least locally). If 
d,pS ^ 0, then the vector field tangent to these trajecto¬ 
ries is given by 




(31) 


The velocity fields Vi and V 2 are generically not the same. 
Using ([^ and ([^, to evaluate the integral in (29), we 
have that 


V 2 =vi-- [ da 
P J — oo 


p 


d^Q 


2{d^SY L 

on C 

{d^ {pd^S) - {pd^s)} 


(32) 

So a wave function ip will lead to the same trajectories 
only if 

on c 

d^Q + {d^ {pd^S) - (pa„5)} = 0. (33) 


Generically the above condition will not hold. For ex¬ 
ample, consider a positive or negative frequency wave 
function. Then the term in curly brackets is zero (be¬ 
cause the density p then equals \ip\^ and hence satisfies 
([^) and the condition is reduced to d^Q = 0, which is 
generically not satisfied. 

On the other hand, if the wave function ip is ei¬ 
ther only left- or right-moving, i.e., ip = iphioi + p) or 
Ip = ipR^a — p), then the condition (33) is automatically 
satisfied. In this case, the trajectories are the same for 
both Bohmian approaches. In addition, the quantum po¬ 
tential is zero, so that the trajectories are classical. 
(Conversely, if Q = 0, then the wave function p is either 
left- or right-moving and the trajectories are the same for 
both Bohmian approaches.) As shown in section [n} the 
classical trajectories are given by a = zLp + c, with c a 
constant. The positive and negative sign respectively cor¬ 
respond to a right- and left-moving wave function. The 
trajectories reach the singularity a = e“ = 0 respectively 
at p ^ —oo and p —> oo. Note, however, that in our sec¬ 
ond Bohmian approach, we have not introduced a lapse 


^ Note that in section|^ we wrote the wave function as ip = ip^ + 
Ip—, whereas in section |III| we wrote it as 'I' = {ip+,ip—). These 
notations are of course equivalent. 


function and hence we do not have a straightforward def¬ 
inition of proper time. So, it is not immediately clear 
whether the singularities are reached in finite proper time 
or not. Nonetheless, since these trajectories are classical, 
it seems natural to introduce a proper time that agrees 
with that of the classical theory. As such, the conclusion 
is again that the singularities are reached in finite proper 
time. 

This dynamics seems less natural than the one given 
earlier. The reason is of course that p plays the role 
of time in this dynamics. However, there seems to be no 
reason to give it that distinguished role. For example, we 
could equally well have taken a to play the role of time. 
This would have resulted in completely different paths. 
Of course, this is only a simplified model. In a more 
serious model, we expect that some of the matter fields 
will represent time on an effective level. But this should 
follow from analyzing a more fundamental Bohmian dy¬ 
namics, similar to (see e.g. PlSillT]), rather than 
be postulated a priori. 


V. SINGULARITIES AND BOHMIAN 
MECHANICS 


In this section, we consider the possibility of singu¬ 
larities (zero scale factor) for trajectories given by the 
Bohmian dynamics of section Jig. In particular, using 
the equilibrium distribution (20) as probability distribu¬ 


tion, we calculate the probability for a trajectory to have 
a singularity, for a given wave function. 

We start by considering the asymptotic behavior of the 
Bohmian trajectories. First, suppose that the wave func¬ 
tion has only right-moving components, i.e., '^{a,p) = 
'ki^(a - p) = - p),pR-{a - p)). Then the 

support of the wave function is localized on a < 0 for 
p ^ —oo and a > 0 for p —> +oo. This follows from the 
fact that 


lim / dap^‘^{a,p) 

M 

= lim / da{\pR^+{a-p)\^ + \pR-{a-p)\ 

0-J.-OO 

= lim f dndpR^+ii^)]^+ \pR-{u)\^) 

0-J.-OO 

= ll'I^itll? = l 

and, similarly, for the other limit we have 

pOO 

dap^^{a, p) 

pOO 

I da {\pR^+{a - p)f + \pR-{a - p)f 


(34) 


lim / 

)->. + 00 Jq 


= lim / 
<P—>-+°o Jq 


= lim / di' I 

= ll^flll? = l- 


‘’R,+ 




( 35 ) 












In the same manner, for a wave function that has only 
left-moving components, i.e., ^'(a, 0) = + cj)) = 

(V'L, + (a -f 4>), ipL,-[ct + </<)), we have that 
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poo 

lim / dap'^^{a,(j)) 

0->—OO Jq 

pOO 

= lim / da (\ipL,+ {a + + li’L-ia + i 

4>^-00 Jq 

poo 

= lim / dv {\i;L,+ i^)\^ + 


= ll'I^L|lt = l 


(36) 


and 


.0 

lim / dap'^^ 

<P^+°oJ-oo 

= lim / da (|•^/'L,+ (a-f -f |V'L.-(a-f </>)|^) 
d-^+ooJ-oo 

= lim [ diy {\tPL,+ H\^ + 
d-^+oo J_oo 


= ll'I'L||f = l 


(37) 


so that the support is now localized on a > 0 for ^ —oo 

and a < 0 for (/) —>■ -foo. 

So in summary, we have that 


'd/R has support on 


has support on 


a<0 for^- 

a>0for^- 
a > 0 ioT (j) 


oo 

► —00 


a < 0 for (j> —>■ oo 


This analysis shows that for any wave function d/ = 
'I'i? + 'I'L) its left- and right-moving components will be¬ 
come asymptotically non-overlapping in a-space. As a 
consequence, asymptotically, Bohmian trajectories will 
be determined either by the or 4'/,. Hence, asymp¬ 
totically, the trajectories are classical, given hy a = (p + c 
or a = —(j) + c, see figure (The same holds for the 
Bohmian dynamics of section [nj in the case that the tra¬ 
jectories case be expressed as a{(j)).) 

This asymptotic classical behavior implies that there 
may be four different types of trajectories. Namely, tra¬ 
jectories that represent a universe that: 

1 . starts with a big bang and keeps expanding forever 

2 . keeps contracting until a big crunch 


a 



FIG. 2: Asymptotic classical behavior of the Bohmian tra¬ 
jectories for an arbitrary wave function. 


3. starts with a big bang and ends in a big crunch 

4. undergoes a bounce, i.e., contracts until it reaches 
a minimal size and then expands again. 

Only trajectories of type 4 are non-singular. As we shall 
see, all four types of trajectories may occur. Though, for 
a given wave function, there can not be trajectories of 
both types 1 and 2. This is an immediate consequence 
of the no-crossing property of the Bohmian dynamics 
(which states that trajectories cannot cross at equal times 
in configuration space), since such trajectories must cross 
each other. 

In the previous section, we showed that wave functions 
that are purely left- or right-moving give rise to classical 
trajectories, which are either of type I or type 2, and 
which are singular. There are also wave functions for 
which there are trajectories of type 3 or 4. 

As a simple example, consider a wave function which 
is even or odd in a. For such a wave function, the ve¬ 
locity field V 2 is odd in a. As such, for every possible 
trajectory a{(j)), also —a{(j)) is a possible trajectory. Be¬ 
cause of the no-crossing property, this implies that no 
trajectory will cross the (p-aods.^ In other words, trajec¬ 
tories that have positive value of a at any given time p, 
must have a > 0 at all times and hence avoid the singu¬ 
larity. These trajectories correspond to a bouncing uni¬ 
verse. They asymptotically start with a = —p + Ci, reach 
a minimal value a^m and then asymptotically evolve ac¬ 
cording to a = p + Cf for p ^ oo. See figurefor an ex¬ 
ample. Conversely, trajectories that asymptotically start 
with a = p + Ci originate from a singularity. Eventually 
they reach a maximal vale Omax and then for p ^ oo 
they move again to the singularity according to the tra¬ 
jectory a = —p + Cf. So these trajectories are big bang 
- big crunch solutions. 


There might also be a trajectory that does not display this 
asymptotic behavior and acts as a bifurcation line between tra¬ 
jectories with different possible asymptotic behavior. 


^ Their might also be the trajectory a{(p) = 0 which acts as a 
bifurcation line between trajectories with different asymptotic 
behavior. 
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FIG. 3: Some trajectories for a wave function that is sym¬ 
metric under a —> —a. The trajectories on the left represent 
universes which start with a big bang and end in a big crunch. 
The trajectories on the right represent universes that bounce. 


Let us now analyze the general case and establish the 
probability to have singular trajectories for a general 
wave function 'F = 'I'a-I-'I'l. Call and Pn^i the prob¬ 
ability for the universe to start respectively with a < 0 
(from the left) and a > 0 (from the right) in the limit 
(j) —)■ — 00 . Similarly, P^j and Prj denote the probabil¬ 
ities for the universe to end up respectively with a < 0 
and a > 0 in the limit (j> —>■ oo. Since and do not 
overlap asymptotically, we have 


PL,i = lim 

(i—>■ —OO 


and 


dap'^{a,<l)) = , 

D 

poo 

PR,i= lim / dap'^{a,(l>) = \\-i>L\\l, 

0->-oo Jq 

f’L,/= lim f dap'^{a,(j)) = \\'i>L\\l, 

0-i-oo 

poo 

Prj = lim / dap'^{a,(j)) = \\'Sr\\1, 

<t>^oo Jq 

PL,i = 1 — Plj = Prj = 1 — Pr,z • 


(38) 


(39) 


We can now find the probability for the universe to 
start or to end in a singularity for a given wave func¬ 
tion d>. First, assume that Pr^^ < Prj. In this case, 
trajectories starting from a > Oas^—>■ —oo can not 
end up moving to a < 0 as ^ oo. If a trajec¬ 
tory did do that, then all the trajectories that started 
on the a < 0 would must also end up with a < 0 as 
(/)—)■ oo because of the no-crossing property. This im¬ 
plies that Prj ^ PR^i and hence, because of (39), that 


PR,i ^ Prj I but this contradicts our assumption. Thus, 
in the case that Pr ^ < Prj-, the probability for a bounce 
Is I^bounce PR,^■ In addition, it implies that the prob¬ 
ability for trajectories to start from a singularity and 
to keep expanding is Pexpanding = Prj - PR,i and the 
probability for trajectories to start and end in a singu¬ 
larity (i.e., trajectories with a big bang and big crunch) 
is Precollapsing — Pl,2 {Pr.j Pr,C) — PR,i- Finally, 


the probability for trajectories coming from a > 0 to 
keep contracting towards the singularity is Pcontracting = 

0. Similarly, for the case that Pr^i > Prj, one can 

show that Pbounce — PrJ, Pcontracting — PlJ PL,i, 
Precollapsing — PL,i and Pexpanding — 0- 

In summary, for an arbitrary state 'k, we have 

I^bounce — .I^recollapsing — Ulin(PR ^5 Plj') ~ Ulin(PR j, Prj') , 
Pexpanding — Cnax(PR ^ PlJ, 9) — ^^^{_PrJ Pr^I, 0) , 
Pcontracting — bnax(PR J Pl,!, 0) — bnax(PR j Prj, 0) ■ 

(40) 

This implies that the probability Psinguiarity = 1-Pbounce 
to run into a singularity satisfies 

2 ^ Psinguiarity ^ 1 ; (^^) 

i.e., the probability to run into a singularity is at least 
one half. Maximum probability is reached when there 
are only left- or right-moving components. In that case, 
every trajectory either runs into a singularity or starts 
from a singularity. On the other hand, for a superpo¬ 
sition of left- or right-moving components, there is al¬ 
ways a non-zero probability for a bounce. For a state 
41 , the maximum probability that can ever be attained 
for a trajectory to be non-singular is 1/2, which happens 
when Pr i = Prj. This happens for the examples of the 
symmetric and anti-symmetric states discussed before. 

The trajectories for the first Bohmian approach, that 
was presented in section [l^ are different from those con¬ 
sidered here. However, the trajectories in the first ap¬ 
proach that are expressible as OL{(f) have the same asymp¬ 
totic behavior as those considered here and hence are 
qualitatively the same. However, there is no natural 
probability distribution over trajectories which allows to 
calculate the probability for trajectories to be singular. 


VI. SINGULARITIES AND CONSISTENT 
HISTORIES 

According to the consistent histories approach the 
probability for a singularity is always one. Coarse¬ 
grained histories are either of type 1 or 2. This was shown 
for two-time histories in lai, with the times being the 
infinite past and future. 

In [T], some of us argued that introducing a third in¬ 
termediate time would not lead to a family of consistent 


The example that was used in [SE] to compare the Bohmian ap¬ 
proach to the consistent histories approach in the context of non- 
relativistic quantum mechanics is actually completely analogous 
to the case of a symmetric state in our cosmological model. In 
Biia, the wave function is composed out of two packets that cross 
each other in time. The wave function is considered to be com¬ 
pletely symmetric so that Bohmian trajectories will never cross 
the symmetry axis, while the histories in the consistent histories 
approach correspond to either left- or right-moving paths. 








histories, unless the state is classical (i.e. the state is a 
left- or right-moving wave function). As such, for non- 
classical states, the consistent histories approach would 
be unable to deal with the properties of the universe at 
intermediate times. This would of course undermine the 
consistent histories approach for this cosmological model. 

However, the argument in [T] is incomplete, mainly be¬ 
cause only one particular family of histories was consid¬ 
ered. That is, the class of propositions (which are repre¬ 
sented by projection operators) was taken to be the same 
as that in HE] for all three times. And for this class, the 
histories are not consistent (unless the intermediate time 
is chosen sufficiently early or late). However, this does 
not exclude other possible families of histories for which 
the histories are consistent. We will show that a suitable 
class of propositions can be chosen such that consistent 
histories can be considered for an arbitrary number of in¬ 
termediate times. With this choice, the probability for 
a singularity is always one. This generalizes the results 
in to an arbitrary number of times. 

To setup the consistent histories framework, we will 
use the triplet As explained in section HI 

(7^2, P 2 ,H) would lead to the same results. 

A coarse-grained history is a sequence of regions 
Ai,..., An in a-space at a sequence of times 4>i^... ^(pn- 
These are obtained by considering an exhaustive set of 
regions {Af^^}, = 1,2,..., of a-space for each time 

4>k, k = l,...,n. For each of these regions, there is a 
projection operator 


P^, = [ Pi (da) (42) 

‘‘■k 

defined in terms of the PVM Pi (da). In terms of these 
projection operators, the class operator for a history h = 
(Ai^,...,A")is 

dh=Pl^ U{Pu.Pu-l)PlPh U{Pn-lAn-2) ■■■ 

*n-l 

■■■ Pl2UiP2,Pi)PiiUiPuM, (43) 

where po is some initial time and U^pkTpk-i) = 
g-iff( 0 fc- 0 fc_i) jg unitary time evolution from time 
pk-i to pk- 

In order to associate probabilities to a family of histo¬ 
ries, the decoherence condition 

«0, for h'^h, (44) 


needs to be satisfied. The probability for a history h is 
then given by 

Ph = \\Chnl- ( 45 ) 


Given a wave function 'k (and given that pi is suf¬ 
ficiently far in the past or pn sufficiently far in the fu¬ 
ture), we can always find a collections of sets {Af^^} such 
that the decoherence condition is satisfied. The reason 
is that the wave function is a superposition of a left- 
moving and a right-moving packet. For each time pk, 
we can choose the regions {Aj^} such that support of 
and n are each approximately within one of the 
Aj)^ (but not necessarily with the same ik)- We then 
have that for each time pk that P^^ '^^{a.pk) ~ 0 for 
all but one A((_, which we denote A)). For A|, we have 
that P^k^L{a,pk) « ^L{a,pk)- We then also have that 

ChL^Lia.po) « 4'i(a,(()„), with = (A^,...,A2). In 
other words, there is a coarse-grained history that ap¬ 
proximately follows the track of the packet 4'/,. Similarly, 
there is a coarse-grained history hji that approximately 
follows the track of the packet 'k^j. 

If one of the pk is such that l and 'kfl; are approxi¬ 
mately non-overlapping at that time (for example by tak¬ 
ing pi sufficiently early or pn sufficiently late, since we 
know that the overlap vanishes in the limit p —)■ ±oo.), 
then 

dhL^ia,po) « '^L{a,pn) , dhR^ia,po) « '^R{a,pn) 

(46) 

and 

Ch'i'ia, (/>o) « 0 , for all Hr, Hr. (47) 


Since 
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(48) 


we have that the decoherence condition (44) is satisfied. 
The only two coarse-grained histories with non-zero prob¬ 
ability are Hr and Hr, with probabilities 

Ph,^\\^L\\\, P,,« 114/^11?. (49) 


In order to consider the question whether there is a 
singularity, we should take pi —)■ —oo and pn —t oo. 
Then (46l-(49) hold and, using a similar notation as in 
the previous section, we have 


Pbounce — Precollapsing ~ 0 i 

Pexpanding = PHl ~ , 

f’contracting = ^’lr„«||4'fl|ll- (50) 

The probability for a bounce is negligibly small. It can 
actually be made arbitrarily small by suitably choosing 
the set {A^^}. So according to the consistent histories 
approach, one can always find a consistent family of his¬ 
tories where the probability for a singularity is one. 


The analysis is completely analogous to that of HE]. The reason 
is that, as already noted in footnote |10| we are considering two 
packets that move across each other, just like in the example of 

HE]. 


Eq. ( |48[ ) follows from the fact that a left-moving packet and 
right-moving packet are orthogonal. Since they have no common 
support at the far past they are clearly orthogonal then. Since 
the inner product is preserved over time, they must be orthogonal 
at all times. 
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VII. CONCLUSION 

We have analyzed a Bohmian approach to the 
Wheeler-DeWitt quantization of the Friedmann- 
Lemaitre-Robertson-Walker model. This Bohmian 
approach agrees with the consistent histories approach 
concerning the probabilities for single-time histories. 
However, it makes different predictions for the probabil¬ 
ity of trajectories or histories to have a singularity. In 
the consistent histories approach, at least for the families 
considered in the present paper, the probability for a 
history to have a singularity is one. On the other hand, 
in the Bohmian approach, for generic wave functions 
(i.e., non-classical wave functions), there is a non-zero 
probability for a trajectory to be non-singular and have 
a bounce. 

So, as was already emphasized in [T], where a different 
Bohmian model was considered to make the compari¬ 
son with the consistent histories approach, the question 
whether or not the Wheeler-DeWitt quantization leads to 
singularities depends very much on the version of quan¬ 
tum theory one adopts. 

In this paper, we have only analyzed the Wheeler- 


DeWitt quantization. It would be interesting to also 
study the Bohmian approach to loop quantization. For 
the consistent histories approach to loop quantization, it 
was recently shown that histories do not have singulari¬ 
ties for generic wave functions [3S]. It is unclear whether 
this is also true for a Bohmian approach. 
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